Novel Spin-texture on the warped Dirac-cone surface states in topological insulators 
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We have investigated the nature of surface states in the Bi2Te3 family of three-dimensional topo- 
logical insulators using first-principles calculations as well as model Hamiltonians. When the surface 
Dirac cone is warped due to Dresselhaus spin-orbit coupling in rhombohedral structures, the spin 
acquires a finite out-of-plane component. We predict a novel in-plane spin-texture of the warped 
surface Dirac cone with spins not perpendicular to the electron momentum. Our k-p model calcula- 
tion reveals that this novel in-plane spin-texture requires high order Dresselhaus spin-orbit coupling 
terms. 



PACS numbers: 

Topological insulators (TI) [T]-[9] realize a novel state 
of quantum matter that is distinguished by topologi- 
cal invariants of bulk band structure rather than spon- 
taneously broken symmetries. Its material realization 
in 2D artificial HgTe-quantum wells [7]-[9] and 3D Bi- 
based compounds p!QUT3] has led to a surge of interest 
in novel topological physics throughout the world-wide 
condensed matter community. The surfaces of 3D topo- 
logical insulators or the interface between two materi- 
als with distinct topological invariants host metallic sur- 
face/interface states. Electrons on the surface of an ideal 
topological insulator have an energy-momentum relation- 
ship in the shape of a Dirac cone. In the presence of such 
a single-Dirac-cone surface or interface band, a number 
of exotic quantum phenomena have been predicted [14]- 
[17] • Recent spin-resolved photoemission measurements 
reveal that the spins of the electrons on the surface Dirac 
cone are locked with their momenta, giving rise to helical 
Dirac fermions without spin degeneracy [3] HJ [101 [18H2Q| . 
This one-to-one locking of the electron spin to the mo- 
mentum comes from a combination of strong spin-orbit 
interaction and inversion symmetry breaking at the sur- 
face. Such spin-texture on the surface Dirac cone leads 
to anti-localization properties and plays a central role for 
exotic quantum phenomena. 

Recent angle-resolved photoemission (ARPES) exper- 
iments have observed the surface Dirac cone in a num- 
ber of topological insulators ^ [TTJ |T8]. In particular, 
Bi2Te3 exhibits a Dirac cone with significant hexagonal 
warpingjU [Til UHl [21], i-e. the shape of Fermi surface 
(FS) or constant energy contour evolves from an ideal 
circle to a hexagon to a snowfiake-like shape with increas- 
ing energy. The hexagonal distortion of the cone does not 
change the Berry phase quantifying its topological invari- 
ant, consistent with its topological order |4l[TT]. Fu [21] 
proposed that Dresselhaus spin-orbit coupling can lead 
to hexagonal warping of the Dirac cone in rhombohedral 
structures. The spins acquire finite out-of-t he-plane com- 
ponents through the orbital channel to conserve the net 



value of the Berry's phase. The resulting finite value of 
the out-of-plane component opens up new possibilities for 
unusual phenomena such as the enhancement of interfer- 
ence patterns around crystal defects and a magnetically 
ordered exotic surface. 

In this Letter, we investigate the nature of surface 
states in the Bi2Te3 family of three-dimensional topolog- 
ical insulators using first-principles calculations as well 
as model Hamiltonians. We reproduce a warped surface 
Dirac cone with finite out-of-plane spin component, and 
further predict a novel in-plane spin-texture with spins 
not perpendicular to the electron momentum. We de- 
velop di k - p model to describe this novel in-plane spin- 
texture which requires high order Dresselhaus spin-orbit 
coupling terms, which has not been previously proposed. 
This is a generic property for surface states in strong 
spin-orbit coupling materials, not limited to topological 
insulators. 

Bi2Te3 is a semiconductor with a rhombohedral crystal 
structure (space group R3m). The unit cell contains five 
atoms, with quintuple layers ordered in Te(l)-Bi-Te(2)- 
Bi-Te(l) sequence. Due to the weak bonds between two 
Te(l) atoms adjacent quintuple layers, the (111) surface 
with Te(l) termination can be easily obtained and is the 
surface usually studied in experiments. To obtain the 
surface states, we use a symmetric 30-atomic-layer slab in 
a hexagonal unit cell. The theoretical band calculations 
are performed with the full-potential linear augmented 
plane wave (LAPW) method using the WIEN2K package 
^22j . The generalized gradient approximation (GGA) of 
Perdew, Burke and Ernzerhof [23] was used to describe 
the exchange-correlation potential. Spin-orbit coupling 
was included as a second variational step using scalar- 
relativistic eigenfunctions as basis. After self-consistent 
charges and potentials are obtained in the symmetric 
slab, we add an extra potential on one side of the surface 
for band computations to remove the degeneracy that 
arises from the two opposite sides of the surface. 

We show the computed surface bands (red) as well as 
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FIG. 1: (Color Online) (a) Electronic structure of Bi2Te3. The surface states are shown as red lines, projected bulk bands by 
shaded area. The inset is the area shaded in gray, (b) Energy dependence of Bi2Te3 surface Dirac cone illustrating hexagonal 
warping. The magnitude of the out-of-plane spin component is indicated by the colors - see color bar. (c) In-plane spin 
direction of the surface Dirac cone. 



the projected bulk band structure (green) in Fig. 1(a) 
along high symmetry directions in the surface Brillouin 
zone (BZ). Note that the lower Dirac cone has a particu- 
larly distorted shape (inset in Fig. 1(a)). The insulating 
phase of the bulk is predicted with a 45 meV indirect gap, 
while the experimental value is 165 meV [4 . This under- 
estimation of band gap is a typical problem of GGA. A 
metallic surface band arises due to the strong spin-orbit 
coupling. Between the f -point and the M-point, an odd 
number of surface bands cross the Fermi level (Ef). This 
is the hallmark of the topological insulator. Topolog- 
ical non-triviality can also be proved by calculating the 
time-reversal invariant Z2 value using wave-function par- 
ity analysis, proposed by Fu and Kane [2 . The value of 
Z2 is found to be -1 due to a single band inversion at the 
T-point. 

The topological surface states form a 2-\-l dimensional 
Dirac-cone-like band dispersion in energy-momentum 
space as shown in Fig. 1(b). The shapes of the energy 
contours for the upper Dirac-cone (projected onto the 
E = plane in Fig. 1(b)) are circles for the low energy 
region, but become hexagonal at ~150 meV above the 
Dirac point. The hexagonal warping becomes more pro- 
nounced at higher energy where a snowflake-like shape is 
obtained around 200 meV. These surface states have been 
observed in ARPES experiments ^ [HI |21] and the 
hexagonal warping effect was described in a k - p model 
by FupT|. Recent Fourier-transform scanning tunneling 
microscope (FT-STM) experiments have shown the non- 
trivial interference pattern due to nesting of the hexago- 
nally warped energy contours [24]. When the Dirac cone 
is warped, it must carry a finite out-of-plane spin com- 
ponent. The out-of-plane spin component has three- fold 



symmetry with up and down spin alternation around a 
circle, whereas the band dispersion displays six-fold sym- 
metry. Since the mirror plane is along the f — M direction 
(chosen as x-direction), the out-of-plane spin component 
has to be zero along this direction as shown in Fig. 1(b). 
The in-plane spin texture in Fig. 1(c) demonstrates a new 
phenomenon, non-orthogonality between spin and mo- 
mentum. The helical nature of the spin is preserved at 
low energies near the Dirac point, where the spins remain 
perpendicular to k. When E is increased, a significant 
departure from this ideal helical behavior is observed, al- 
though the spin and momentum remain orthogonal to 
each other along the high symmetry directions F-M and 
F-K. Similar features were also observed in recent pho- 
toemission experiments on Bi2Se3 [25] and it was noted 
that the existing k ■ p theory could not explain the non- 
orthogonality between electron spin and momentum. To 
understand the origin of this non-trivial behavior of the 
spins, we study the surface states of Bi2Te3 using a /c • p 
model with higher order corrections in the following para- 
graph. 

In our formalism, the effective Hamiltonian for the sur- 
face bands is expanded up to fifth order in k, taking the 
following form, 

H{k) = Eo{k)^Vk{k^ay-kya^)^^{kl^k^_)a, 

+ <(4a+-fe^a_), (1) 

where ^o(^) = A^^/2m*, k± = kx ^ iky^ a± = ax ^ i(Jy 
and the cr^ are Pauli matrices. The form of H(k) is suit- 
able for describing the [111] surface band structure near 
the F point in the surface Brillouin zone of Bi2Te3 fam- 
ily TIs and it is invariant under time-reversal and C^y 
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symmetries. Eo^k) generates particle-hole asymmetry 
and the term, HQ{k) = Vk{kx(Jy — kya^), describes an 
isotropic 2D helical Dirac fermion. The Dirac velocity, 
Vk = '^(1 + (^k^ + Pk"^)^ where v is the Fermi velocity, 
contains a fourth order correction term. More impor- 
tantly, the /c^-term, = ^(^+ + k^-)crz, leads to the 
hexagonal warping of the Fermi surface. The hexago- 
nal warping parameter also has a second order correction 
term, A/e = A(l + 7/^^). The last term in Eq. [l] describes 
fifth-order spin-orbit coupling at the surface of rhombo- 
hedral crystal systems, and, to the best of our knowledge, 
it has not been reported before. The purpose of intro- 
ducing this term is to explain the non-orthogonality of 
the spins, which cannot be reproduced by the /c^-term. 
We determine the parameters by fitting the model spin 
directions to the first principles results. The resulting 
values of the parameters are recorded in table |lj 

TABLE I: Parameters used \n k ■ p calculations. 



Parameters 


BisTes 


m* (in eV"' • A"^) 
V (in eV • A) 
a (in eV • A^) 
A (in eV • A^) 
13 (in eV • A^) 
7 (in eV • A^) 
C (in eV • A^) 


35.21 
-0.0005 
14.82 
-0.04 
-1.26 
-0.001 
0.35 



The in-plane spin polarization obtained from our 
model calculations agrees with the first-principles result 
by capturing the non-orthogonality between electron spin 
and momentum. For a quantitative analysis, in Fig. |2] 
(a) we plot the angle {0) between the spin direction ob- 
tained from k • p theory and kx direction as a function 
of the azimuthal angle, (j) at an energy —30 meV. For 
ideally helical spins, 9 would follow the blue dashed line, 
but our calculated spins follows the red curve displaying 
a complicated pattern of spin canting. The calculated 
angle is equal to the ideal angle only at high symmetry 
directions where the two curves intersect. In Fig. [2] (b) 
we highlight the deviation of the angle between spin and 
momentum from orthogonality by plotting the angle of 
deviation, ^ as a function of (j). 5 (defined in the inset) 
oscillates about as we move along a constant energy 
contour in the higher /c-region, where the /c^-term dom- 
inates making non-orthogonal feature more prominent. 
Here also we find that 5 is exactly equal to along F- 
M and F-K directions, indicating that electron spin and 
momentum remain orthogonal along high symmetry di- 
rections in /c-space. 

Figure 3 shows the surface state dispersions obtained 
from the k • p model along the F-M (b) and F-K direc- 
tions (d). The dispersion of the surface state starts out 
linear in momentum /c, as expected for a massless Dirac 
cone, but both lower and upper branches of the disper- 
sion curve rapidly deviate from linearity, becoming con- 
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FIG. 2: Quantitative analysis of non-orthogonality of electron 
spin and momentum on the surface Dirac cone of Bi2Te3. (a) 
Angle between spin s, obtained from k • p theory and kx is 
plotted against azimuthal angle, cj) at an energy —30 meV. The 
blue dashed curve and solid red curve represent ideal and non- 
ideal helical behavior, respectively, (b) Angle of deviation, 5 
(inset) vs. (j) plot. 



vex functions of /c, similar to the GGA results shown 
in Fig. 3 (a) and (c). Both STM and ARPES studies 
of Bi2Se3 capture similar behavior of the surface states 

In conclusion, we have shown that the spin of an elec- 
tron on the surface of a strongly spin-orbit coupled ma- 
terial does not have to be orthogonal to the electron mo- 
mentum due to higher order Dresselhaus spin-orbit cou- 
pling. This generic property should exist in both topo- 
logical insulators and non-topological materials (metals 
or insulators). For example, pure bismuth is topologi- 
cally trivial and may have novel spin-texture. The non- 
orthogonality of spin and momentum should be observ- 
able on the warped surface Dirac cone in the topological 
insulator Bi2Te3 where the spin-texture is crucial for a 
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FIG. 3: Dispersion of the surface states of Bi2Te3 from (a),(c) 
first-principles calculations and (b),(d) k ■ p model. 



number of exotic quantum phenomena. 
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